CLNS 00/1708 
CLEO 00-27 



Search for CP violation in D° -> K^tt , D° -> ttV and D° -> KgKg 

decays 

CLEO Collaboration 
(February 7, 2008) 

Abstract 



We have searched for CP-violating asymmetries in neutral charm meson de- 
cays in 13.7 ftr 1 of e + e~ collision data at^« 10.6 GeV with the CLEO 
detector. The measured asymmetries in the rate of D° and D° decays to 
K^tt , vr 7r° and KgKg final states are (+0.1 ± 1.3)%, (+0.1 + 4.8)% and 
(-23 ± 19)%, respectively. 
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Measurable CP violating phenomena in strange |I].|2[] and beauty [§-|5| mesons are the 
impetus for numerous current and future experiments |6|-f| that are expected to challenge 
the standard model (SM) description of the weak interaction. In contrast the SM predictions 



for CP violation in the charm meson system of (9(0.1%) |T0[ are probably not attainable 



by current experiments, although recent conjectures [O] indicate that direct CP violating 
effects may be as large as (9(1%). Thus an observation of CP violation in charm decays 
exceeding the percent level would be strong evidence for non-SM processes. 



Previous searches for mixing- induced [12| or direct CP violation |T3|JT4| in the neu- 
tral charm meson system have set limits of ~ 30% or a few percent, respectively. We 
present results of searches for direct CP violation in neutral charm meson decays to pairs 
of light pseudoscalar mesons: Kg7r°, 7r°7r° and KgKg. Decays to the latter two final states 
are Cabibbo-suppressed, thus enhancing the possibility that interference with non-SM am- 



plitudes could produce direct CP violation. A previous search by CLEO ||T3[ for direct 
CP violation in D° and D° decays to KgV in 2.7 fb™ 1 of e + e~ collision data established 
A(KI-k°) = (-1.8 ± 3.0)% with the definition 

^ n = r(D°^/)-r(DQ^/) 

r( D o ^ /) + r(D° -►/) K J 

where / is the final state. 

The current results are based upon 13.7 fb -1 of e + e~ collision data at ^/s ~ 10.6 GeV 
accumulated with two configurations of the CLEO experiment at the Cornell Electron Stor- 
age Ring (CESR). Approximately one-third of the data were accumulated with the CLEO 
II configuration |TJj that consists of 3 nested cylindrical wire chambers surrounded by a 
CsI(Tl) electromagnetic calorimeter immersed in a 1.5T solenoidal magnetic field. The 3.5 
cm radius beam pipe and innermost wire chamber were replaced by a 2 cm radius beam pipe 
and a three-layer double-sided silicon vertex detector in CLEO II. V ||16|| . In addition the gas 



mixture in the main drift chamber was changed from an argon:ethane to helium:propane [17| 
mixture for improved charged particle momentum resolution and efficiency 

The Monte Carlo simulation of the CLEO II and CLEO II. V detector configurations was 
based upon GEANT . Simulated events were processed in the same manner as the data. 
We used a simulated sample of e + e~ — > qq (q = u, c, s, d) events representing a luminosity 
comparable to that of the data to determine selection criteria and investigate some systematic 
effects. Systematic uncertainties in the asymmetry measurements are determined from the 
data when possible. 

The charge of the slow pion produced in the decay D* + — > D ^^ identifies the flavor 
of the neutral charm meson at production (charge conjugation is implied throughout unless 
explicitly stated otherwise). Candidate 7r^" ow must be well-reconstructed tracks originating 
from a cylinder of radius 3 mm and half-length 5 cm centered on the e + e~ interaction point. 
A minimum momentum requirement on D° candidates of 2 GeV/c sets a lower limit on the 
7r^ ow momentum of approximately 95 MeV/c. 

Candidates for the decay Kg — ■> 7r + 7r - are formed from opposite-sign pairs of charged 
particles within 8 (9) MeV of the known K s mass @ for the K s Kg (K s tt ) final state. The 
reconstructed K s decay vertex must be separated from the interaction point by at least 3 
standard deviations (a) where a is calculated from the track covariance matrices. In addition 
the x 2 f° r each K s daughter track to originate from the interaction point is required to be 
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larger than 2.5. The latter two requirements not only suppress combinatorial background 
arising from random combinations of tt + tc~ candidates but also diminishes the contribution 
of D° -> h + h-K° s and D° -> h+hr^ (h = n,K) backgrounds to the K° S K° S and K° s ir° final 
states, respectively. 

Neutral pion candidates are formed from pairs of electromagnetic showers in the CsI(Tl) 
calorimeter unassociated with a charged track. Showers in the barrel (end cap) region of the 
calorimeter must exceed 30 (50) MeV to be considered as a 7r° daughter candidate where the 
barrel is the region of the calorimeter at least 45° from the e + e~ collision axis. The invariant 
mass of 7T° candidates must lie within 20 (18) MeV of the known ir° mass [14| for the 7r°7r° 
(KgV ) final state. 

Neutral pion and Kg candidates are kinematically constrained to the 7r° and Kg mass [14 



and combined to form D° candidates. The mass constraints on the D° daughter candidates 
improves the D° mass resolution by 8%, 5% and 17% for the KgV , 7r°7r° and KgKg final state, 
respectively. A final requirement is placed on cos 9d where 9d is the angle in the D° rest frame 
between the tt° (Kg) direction and the D° flight direction for 7r°7r° (K s 7r° and K S K S ) decays. 
Combinatorial background due to low momentum ir° and Kg candidates is peaked towards 
[cos 6^1 = 1 and the two-body decays of the spinless D° have a flat distribution in cos^. 
We require cos6 d to be in the range [-1.00, +0.95], [-0.875, +0.875] and [-0.96, +0.96] for 
K g 7T°, 7T 7r° and K g K s final states, respectively. 

D° candidates are selected by requiring M, the reconstructed D° candidate mass, to be 
within 50, 65, and 18 MeV of the known D° mass |TJ for K s 7r°, 7r°7r° and K S K S final states, 



respectively. The Q distributions of the candidates in the three decay modes are shown in 
Figures [I], ^ and [| respectively, where Q is the energy release, Q = M(D°7r^ ow ) — M — M T +, 
M(D 7r^ ow ) is the D°7r^ ow invariant mass and M n + is the charged pion mass |L4|]. A prominent 
peak indicative of D* + — > D°7r,^ ow decays in observed in all three distributions. 

The sum S of the number of D° and D° candidates to a given final state (the denominator 
in Eqn (fj) ) is determined by fitting the background in the Q distribution. The background 
shape is approximated as a non-relativistic threshold function with first and second order 
relativistic corrections B(Q) = aQ 1 / 2 +bQ 3 / 2 + cQ 5 ^ 2 and the signal region Q e [3.3, 8.3] MeV 
is excluded from the fit. The interpolated background in the signal region is determined from 
the fit and subtracted from the total number of D° and D° candidates to determine S. For the 
three decay modes under investigation, we obtain iS(K s 7r°) = 9099 + 151, iS(7r°7r°) = 810 + 89 
and iS(K s K s ) = 65 + 14 where the quoted uncertainty includes the uncertainty due to the 
background interpolation. The numerator in Eqn. [TJ is determined from the difference in the 
number of D° and D° candidates in the region Q e [3.3, 8.3] MeV. The measured raw asym- 
metries are (+0.0 ± 1.1)%, (+0.1 ± 4.8)% and (-14 ± 14)% for the KgV°, ttV and K S K S 
final states, respectively. This method of determining the asymmetry implicitly assumes 
that the background is symmetric. As shown in Figures |l|, |2| and |3|, the Q distributions are 
indeed symmetric outside the region Q 6 [3.3, 8.3] MeV to the statistical precision available. 

We have measured the momentum-dependent detector- or reconstruction-induced slow 
pion asymmetry by selecting charged pions from Kg decays using the same selection criteria 
used to select D* + daughters. Since the inner detector material differs for the two configura- 
tions, we take the measured asymmetry as a function of momentum for each configuration 
and weight it by the vr^ ow spectrum from D* + decays. The asymmetry measured for CLEO 
II and CLEO II. V is (-0.20 ± 0.34)% and (+0.18 ± 0.23)%, respectively, where the uncer- 
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FIG. 1. (a) Fitted Q distribution for D° — > Kg7r°. The points with error bars are the data, the 
(barely visible) solid line represents the fitted background and the dashed line shows the interpo- 
lation into the Q signal region, (b) The Q distributions for D° — > Kg7r° (points) and D° — > Kg7r° 
(histogram) candidates. 
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FIG. 2. (a) Fitted Q distribution for D° — > 7r°7r°. The points with error bars are the data, 
the solid line represents the fitted background and the dashed line shows the interpolation into 
the Q signal region, (b) The Q distributions for D — * (points) and D° — > 7r°vr (histogram) 
candidates. 
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FIG. 3. (a) Fitted Q distribution for D° 
the solid line represents the fitted background and the dashed line shows the interpolation into the 
Q signal region, (b) The Q distributions for D° K S K S (points) and D° -» K§K§ (histogram) 
candidates. 



KgKg. The points with error bars are the data, 
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tainties are statistical only. The overall asymmetry, weighted by the luminosity accumulated 
with each configuration, is (+0.05 ± 0.19)%. We apply a correction to the measured raw 
asymmetries for D° decays for this small bias. 

We have also investigated the possibility of a bias in the fitting method by repeating 
the procedure on samples of simulated events. Q distributions of D° and D° candidates 
with known asymmetries ranging between ±50% and with statistics comparable to that of 
the data were constructed from simulated events. In addition we measured the asymmetry 
by fitting the Q distributions to obtain the number of D° and D° candidates separately, 
instead of the procedure described above. We also used an alternative parameterization 
of the background that was found to accurately model the background shape in simulated 
events: B(Q) = e c ®Q n where c and n are fitted parameters. The magnitude of the average 
bias observed in these samples was < 0.1% and we estimate the systematic uncertainty from 
the fitting procedure to be 0.5% (approximately twice the RMS of the asymmetries measured 
in the simulated samples). 

Finally we must take into account the possible effects of asymmetries in other D° decay 
modes that have mistakenly been reconstructed in one of the three selected decay modes. 
We can write the measured asymmetry A m as A m = (A s + rAb)/{l + r) where A s is the 
asymmetry of the selected D° decay mode ("signal"), Ab is the asymmetry of other D° 
modes ("background") that are reconstructed as the signal mode and r is the ratio of the 
rates of the background to the signal. Background from decays of the form D° — > Kg7r°7r°, 
7r + 7T~7r 7r°, 7r 7r 7r° or K s K s 7r can not contribute because their reconstructed mass is ap- 
proximately one pion mass below the D° mass and well outside the allowed range for M. 
Background from D° — > h + h~7r° (D° — > h + h~K^ and D° — > h + h~h + h~) can contribute 
to the reconstructed Kg7r° (K s K g ) yield when the h + h~ mass is near the Kg mass. There 
are no backgrounds of this sort that can contribute to the 7r°7r° yield. The magnitude and 
asymmetry of this background can be directly measured using K s — > 7r + 7r~ candidates with 
M(it + it~) in sidebands either just below or above the M(7r + 7r~) mass range for standard Kg 
candidate selection. Using the same analysis procedure for these sideband candidates, we 
measure r(K s K s ) = 0.16±0.11 and A(K S K S ) = (+40±42)% and r(K°7r°) = 0.03±0.02 and 
^4fc(Kg7r°) = (—5.5 ± 5.1)% where the uncertainties are statistical only. The relative back- 
ground rate for KgKg is substantially higher than that for Kg7r° because the primary decay 
mode contributing to the KgKg background is the Cabibbo-favored D° — > K s 7r + 7r~ while 
the main contributors to K g 7r° are the Cabibbo-suppressed D° — > 7r + 7r~7r° and kinematically 
asymmetric D° — > K~-k + tt decays. 

Correcting the measured raw asymmetries for the slow pion reconstruction bias and the 
rate and asymmetry of the background, we obtain ^4(K s 7r°) = (+0.1 ± 1.3)%, ^4(7r 7r ) = 
(+0.1 ±4.8)% and *4.(K s Kg) = (—23 ± 19)% where the uncertainties contain the combined 
statistical and systematic uncertainties. All systematic uncertainties, except for that assigned 
for possible bias in the fitting method, are determined from the data and would be reduced in 
future higher luminosity samples. All measured asymmetries are consistent with zero and no 
indication of significant CP violation is observed. The former measurement is a substantial 
improvement over the previous CLEO measurement |H| and supersedes it. This is the first 
measurement of the asymmetry in D° decays to the Cabibbo-suppressed final states K s Kg 
and 7r 7r°. 
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